Exciting Collective Oscillations in a Trapped ID Gas 
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We report on the realization of a trapped one dimensional Bose gas and its characterization by 
means of measuring its lowest lying collective excitations. The quantum degenerate Bose gas is 
prepared in a 2D optical lattice and we find the ratio of the frequencies of the lowest compressional 
(breathing) mode and the dipole mode to be (lub/lod) 2 — 3.1, in accordance with the Lieb-Liniger 
and mean-field theory. For a thermal gas we measure {lub/^d) 2 — 4. By heating the quantum 
degenerate gas we have studied the transition between the two regimes. For the lowest number of 
particles attainable in the experiment the kinetic energy of the system is similar to the interaction 
energy and we enter the strongly interacting regime. 

PACS numbers: 05.30.Jp, 03.75.Kk, 03.75.Lm 



o 

CO 



T3 
S3 
O 

o 



> ■ 
■ 

o : 
\o . 
t> ■ 
o ■ 
co : 
o . 

-i— > ■ 
a ■ 



i 

-a 

o 
o 



X 



An ultra cold Bose gas in one spatial dimension is dif- 
ferent from its two or three dimensional counterparts. 
One striking example is that Bose-Einstein condensation 
does not occur at finite temperature in a homogeneous 
one dimensional system. In an interacting Bose gas the 
constraint to one dimension leads to another remarkable 
and counterintuitive property. With decreasing atomic 
density the interactions become increasingly dominant 
and the character of the system changes. Assuming delta- 
functional interactions, exact solutions have been found 
for the ground state and the excitation spectrum of a 
homogeneous one dimensional Bose gas jjj, |2( • Sparked 
by the prospect that this unique model in many-body 
quantum physics could become experimentally accessi- 
ble there has recently been a wave of theoretical interest 
in trapped ID gases. Assuming elongated trapping ge- 
ometries in which the radial atomic motion is confined to 
zero point oscillation different physical regimes could be 
identified 0, S IE IE S S @ • The stringent requirements 
for testing these models have so far not been reached in 
experiments. 

A trapped ID gas is characterized by a single param- 
eter 7 which is the ratio between the interaction energy 
and the kinetic energy of the ground state: 7 = p^j; > 
with m being the atomic mass, gm the ID coupling con- 
stant and niD the ID density. For high densities (7 -C 1) 
the system is weakly interacting, a regime which can 
be described by mean-field theory, and in harmonically 
confined ID systems Bose-Einstein condensation is pos- 
sible. When the ID density is lowered the kinetic en- 
ergy of the ground state is reduced and may get smaller 
than the interaction energy thereby transforming the gas 
into a strongly interacting system, where the longitu- 
dinal motion of the particles is highly correlated. For 
7 > 1 the correlation length l c = %/ ^JmrimgiD becomes 
shorter than the mean interparticle distance. In the limit 
7 — ► 00, often referred to as the Tonks-Girardeau regime, 
a Bose gas acquires fermionic properties due to the strong 
repulsive interactions. The different regimes associated 
with the parameter 7 are characterized by the excita- 



tion spectrum which can be probed by measuring the 
frequencies of collective oscillations |4j, |3, |<J . In the two 
limiting cases of weak and strong interactions the ground 
state properties and the excitation spectrum are explic- 
itly known, whereas in the crossover regime many un- 
known issues remain to be resolved for the confined sys- 
tem. 
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FIG. 1: The geometry and size of trapped ID gases in a two 
dimensional optical lattice. The spacing between the ID tubes 
in the horizontal and vertical direction is 413 nm. 



The ID regime is reached, when the condition 
fj,, ksT <C fojJ r 



(1) 



is fulfilled, where u r denotes the radial trapping fre- 
quency and thus the strength of the radial harmonic 
confinement, \i the chemical potential, and T the tem- 
perature. If the ground state extension in the radial 
direction a r = \f%] muj r is much larger than the char- 
acteristic radius of the interatomic potential, the ID 
coupling strength gm can be expressed in terms of the 
3D scattering length a through guj = — with 
aio — —a 2 , /a being the ID scattering length 3]. The 
necessity for low densities and therefore low atom num- 
bers makes the experimental quest for ID gases challeng- 
ing. A simple estimate using the Lieb-Liniger chemical 
potential /1 = 1hw r n\r)a shows that for a degenerate ID 
quantum gas in the mean field regime the density must 
obey riio <C 1/a and the number of particles must be 
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N <C L/a ~ 10 3 — 10 4 , where L is the characteristic 
length of the system, e.g. the Thomas-Fermi radius. 

There has been significant progress towards the real- 
ization of trapped ID atom ic g ases over the pas t years. 
Both, in a 6 Li/ 7 Li mixture and in 23 Na [ll|, quan- 
tum degenerate gases have been created in very elongated 
traps and features of one dimensional condensate expan- 
sion were observed. Considering only the condensed frac- 
tion a chemical potential of /i > 0.5 Tiuo r was attained 
in these experiments. However, the thermal component, 
which was a substantial portion of the gas, was in a 3D 
configuration (ksT > huj r ) leaving the whole sample 
in an interesting crossover regime. In a similar exper- 
imental regime 7 Li Bose-Einstein condensates with at- 
tractive interparticle interactions were launched into ID 
matter waveguides forming bright matter wave solitons 
[l^ . Moreover, a Bose-Einstein condensate of rubidium 
has been loaded into the ground state of a two dimen- 
sional optical lattice, where the transverse oscillations 
have been frozen out but the tunneling rate between the 
tubes exceeded the axial trapping frequency, resulting in 
an array of strongly coupled tubes • 

In our experiment we have realized both quantum de- 
generate and thermal one dimensional atomic gases with 
the condition (1) being well fulfilled: for all our experi- 
ments ksT/haJr < 6x 10~ 3 and fi/~huj r < 0.1. The atoms 
are prepared in a 2D optical lattice which offers the ad- 
vantage of an extremely tight radial confinement of only a 
fraction of the optical lattice wavelength. Moreover, the 
geometry (see Fig. QJ makes it possible to study many 
copies of the ID system at the same time, thereby cir- 
cumventing problems arising from the detection of very 
low particle numbers. The parameter 7 ranges approxi- 
mately from 0.4 to 1, which is at the crossover from the 
mean field to the strongly correlated regime. 

In the experiment, we collect up to 2 x 10 9 87 Rb 
atoms in a vapor cell magneto-optical trap. After po- 
larization gradient cooling and optical pumping into the 
\F = 2, nip = 2) hyperfine ground state the atoms are 
captured in a magnetic quadrupole trap. Subsequently, 
we magnetically transport the trapped atoms u sing a se- 
ries of partially overlapping quadrupole coils |lj] over 
a distance of 40 cm into an optical quality quartz cell, 
which is pumped to a pressure below 2 x 10 _11 mbar. 
The transfer time is 1.5 s and we have observed transfer 
efficiencies larger than 80% with no detectable heating 
of the cloud. Finally, the linear quadrupole potential 
is converted into the harmonic and elongated potential 
of a QUIC-trap 0. Radio frequency induced evapo- 
ration of the cloud is performed over a period of 25 s, 
during which the confinement is adiabatically relaxed to 
■jj±_ = 2n x 120 Hz in the radial and uu = 2tt x 20 Hz 
in the axial direction. We produce almost pure Bose- 
Einstein condensates of up to 3 x 10 5 atoms. After con- 
densation we adiabatically change the trapping geometry 
to an approximately spherical symmetry with trapping 



frequencies of ui x — 2ir x 17 Hz, uj y — 2ir x 20 Hz, and 
uj z = 2ir x 22 Hz. This reduces the peak density by a 
factor 4 and allows us to load the optical lattice more 
uniformly. 

The optical lattice is formed by two retro-reflected 
laser beams, which are derived from laser diodes at a 
wavelength of A = 826 nm. At the position of the con- 
densate, the laser beams overlap perpendicularly with or- 
thogonal polarizations and are focused to a circular waist 
(l/e 2 -radius) of 105 /im and 120 /im, respectively. The 
frequencies of the two beams are offset with respect to 
each other by 152 MHz. The optical potential depth U 
is proportional to the laser intensity and is conveniently 
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expressed in terms of the recoil energy E rec = with 
k = =?■. The two dimensional optical lattice produces an 
array of tubes, tightly confined in the radial direction and 
spaced by the periodicity of the lattice d = A/2 (see Fig. 
HJ. The radial confinement of the tube is characterized 
by the radial trapping frequency hu) r ~ 2^/sE rec , where 
s = U I ' E rec expresses the potential depth in units of E rec . 
The axial confinement of the tubes is determined by the 
waists w and intensities of the laser beams and character- 
ized by the trapping frequency %w z ~ 2y/s^E rec . We 
have calibrated the potential depth of each of the op- 
tical lattice laser beams by measuring the frequency of 
small amplitude dipole oscillations along the axis of the 
laser beam. From the oscillation frequency we deduce 
the effective mass m* /m at the quasi- momentum q = 
in the band structure, which is a measure of the potential 
depth of the optical lattice 0, 0] . The calibration error 
is estimated to be < 10%. 

Adiabatic loading into the ground state of the optical 
lattice was achieved by ramping up the laser intensity to 
U = 30 E rec with an exponential ramp using a time con- 
stant t = 75 ms and a duration of to = 150 ms. We have 
verified experimentally that all atoms are loaded into the 
lowest Bloch band of the optical lattice by ramping down 
the intensity of the lattice laser beams adiabatically and 
observing in the time-of-flight image, after release from 
the magnetic trap, that only the lowest Brillouin zone 
was occupied [lflj . 

The ID systems in the optical lattice are not per- 
fectly isolated, but the tubes are coupled by the tun- 
nelling matrix element J. For sufficiently deep lattice po- 
tentials, the tunnelling becomes exponentially small and 
contributes only a small correction of order J//i to the ID 
characteristics in the individual tubes. If J/fi <C 1 locally 
the gas acquires ID properties and can be well described 
by a local Lieb-Liniger model, even though the whole 
sample is three dimensional 0. Experimentally we ob- 
serve the disappearance of the matter wave interference 
pattern with increasing lattice depth when the atoms are 
suddenly released from the optical lattice. Higher order 
momentum peaks (±2hk, zLATik, ...) are usually observed 
at lower laser intensities We attribute this loss of 

coherence between the individual tubes to the very small 
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tunnel coupling at large lattice depths, which is too small 
to stabilize the global phase coherence. 

The different regimes of the ID gas in the optical lattice 
can be determined from the excitation spectrum. The 
frequency ratio between the lowest compressional mode 
(breathing mode) and the dipole oscillation (lub/ud) 2 is 
a sensitive measure, both for isolated ID systems |7j and 
for atoms in an optical lattice For a degenerate gas 
in the ID mean field regime, one expects (ojb/i^d) 2 = 3, 
whereas in the Tonks-Girardeau regime (ub/ijJd) 2 = 4. 
The latter frequency ratio is the same for both a thermal 
gas and a gas of degenerate, noninteracting fermions. In 
contrast, for a three dimensional elongated condensate in 
the mean field regime the ratio of the oscillation frequen- 
cies is (ujb/uJd) 2 = 5/2 HIIU El. 



ballistic expansion by absorption imaging. The density 
distribution of the atoms was fitted by a Gaussian to ex- 
tract the position and width of the cloud. Figure 2 shows 
a data set of a dipole oscillation (Fig. 2a) and a breath- 
ing mode (Fig. 2b). In order to extract the frequencies of 
the modes, we have fitted an exponentially decaying sine 
function to the position of the cloud and an exponentially 
decaying sine function plus a linearly increasing term to 
the rms axial width of cloud. The latter accounts for the 
observation that there is a slight increase in the width 
of cloud with longer hold times, possibly due to techni- 
cal noise. The damping coefficients varied between 0s _1 
and 40 s^ 1 for the dipole oscillations and between 3s _1 
and 60s -1 for the breathing mode. We found that the 
damping coefficients depend critically on the alignment 
of the optical lattice. 
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FIG. 2: Dipole oscillation (a) and breathing mode (b) of a 
quantum degenerate one dimensional Bose gas. For this data 
set an almost pure Bose- Einstein condensate with N = (9.8 ± 
0.8) x 10 4 atoms was loaded into the optical lattice and imaged 
after 15 ms of ballistic expansion. From the fits we obtain 
lu d = 2tt x (84.6 ± 0.4) Hz and lu b = 2tt x (152.6 ± 2.0) Hz. 

We have measured the frequency of the collective exci- 
tations of the atoms in the optical lattice. The breathing 
mode was excited by sinusoidal intensity modulation of 
the optical lattice with an amplitude of 4:E rec for five 
cycles and a frequency of 150 Hz, which is close to but 
not matching the expected frequency of the breathing 
mode. At the end of the modulation period a short (1 
ms) magnetic field gradient was applied along the sym- 
metry axis of the ID tubes to induce a dipole oscillation 
of the condensate in the axial trapping potential. Af- 
ter a variable evolution time in the combined optical and 
magnetic trapping potential all confining forces were sud- 
denly switched off |22j and we detected the atoms after 
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FIG. 3: The measured frequency ratio (lob/^d) 2 for a Bose 
condensed ID gas (solid symbols). The solid curve is the 
theoretical prediction from [Sj using the measured dipole fre- 
quency along the ID tubes ui z — 2ir x 84 Hz and the frequency 
u> of the slowly v arying confining potential in the transverse 
direction with ^/m/m*uj — 2ir x 4 Hz. Averaging all mea- 
surements, independent of the atom number, we obtain a fre- 
quency ratio of (ojb/ojd) 2 = 3.15 ± 0.22. For a ID gas of 
thermal atoms we find {uib/^d) 2 = 4.10 ± 0.08 (open cir- 
cles). The depth of the optical lattice is 30-E rec . The error 
bars reflect only the statistical uncertainties on the total atom 
number and fit errors on the frequencies. 

Figure 3 shows the measured ratio (lob/^>d) 2 for a pure 
3D-Bose-Einstein condensate loaded into a 2D optical 
lattice for various total atom numbers. For all data points 
there was no discernible thermal cloud which allows us to 
estimate for the temperature T/T c ^d < 0.3, where T c ,3zj 
denotes the critical temperature for Bose-Einstein con- 
densation in the final magnetic trapping configuration. 
We compare the measured ratio (lvb/^d) 2 to the theo- 
retical prediction of (solid line) and find good agree- 
ment over the wide range of atom numbers investigated. 
For the lowest total atom numbers N = 1.7 x 10 4 the pa- 
rameter 7 reaches unity in the central tube of the lattice, 
indicating that we are in the crossover region from the 
ID mean field regime to the Tonks-Girardeau regime. We 



4 



have estimated the number of atoms in the central tube 
to be 30, assuming that the overall 3D density profile is 
Thomas-Fermi-like and using an effective coupling con- 
stant g which is modified by the optical lattice |l7J, LUt • 
We have also loaded thermal gases into the optical lat- 
tice and obtained an average value (ujb/^d) 2 = 4.10 ± 
0.08 without significant dependence on total atom num- 
ber iV and temperature T over the range of 6 x 10 4 < 
N < 1.6 X 10 5 and 54 nK < T < 91 nK. For thermal 
clouds we have observed that the frequency of the dipole 
oscillations is up to 5% smaller than for the Bose con- 
densed clouds. We attribute this to the larger size of the 
thermal clouds which therefore might experience anhar- 
monic parts of the optical potential. 
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FIG. 4: The measured frequency ratio (ojb/ujd) 2 for a Bose 
gas with 7 x 10 4 atoms after heating in the optical lattice. 
The rms axial width is given after 15 ms of time-of- flight. 
The inset shows the evolution of the axial width vs. hold 
time in the optical lattice prior to excitation of the collective 
modes. 

To study the transition from the ID quantum de- 
generate gas to the ID thermal gas we have prepared 
atomic clouds in the optical lattice with increasing non- 
condensed fraction but constant atom number. We 
have loaded an initially pure Bose-Einstein condensate 
of 7 x 10 4 atoms into an optical lattice and trapped it 
there for different hold times before exciting the collec- 
tive modes. During the hold period the condensate is 
subjected to heating by off-resonant photon scattering 
and possibly technical noise on the trapping fields. We 
find that together with the hold time the axial width of 
the atomic cloud increases (see inset of Fig. 4), whereas 
the radial size is unaffected. Since the timescale for ther- 
malization of the ID gas is unknown it remains unclear 
whether the cloud is in thermal equilibrium and we re- 
frain from calculating a temperature from the rms axial 
width of the cloud. Figure 4 shows the measured ra- 
tio (ujb /^d) 2 as a function of the rms axial width of 
the cloud. For increasing width the ratio (ujb/lud) a P~ 
proaches 4, which is the value for a classical noninteract- 
ing gas. From a simple estimate we deduce that the ID 



gases are in a collisional regime along the axial direction, 
since Nm • (1 — T) ^> 1, where T is the transmission coef- 
ficient for a ID collision of two bosons Q and N\d ~ 70 
is the number of atoms in a ID tube. 

In conclusion, we have realized both thermal and quan- 
tum degenerate gases in one dimension and investigated 
their physics by measuring the low-lying collective exci- 
tations. Our measurements show that the properties of 
the ID ground state are extremely sensitive to thermal 
excitations and finite temperature effects must be taken 
into account when studying ID gases, in particular for 
the identification of the Tonks-Girardeau regime of im- 
penetrable bosons. 
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